Thin-Film Magnetization Dynamics on the Surface of a Topological Insulator 



Yaroslav Tserkovnyak 1 and Daniel Loss 2 

1 Department of Physics and Astronomy, University of California, Los Angeles, California 90095, USA 
2 Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland 

We theoretically study the magnetization dynamics of a thin ferromagnetic film exchange-coupled 
with a surface of a strong three-dimensional topological insulator. We focus on the role of electronic 
zero modes imprinted by domain walls (DW's) or other topological textures in the magnetic film. 
Thermodynamically reciprocal hydrodynamic equations of motion are derived for the DW respond- 
ing to electronic spin torques, on the one hand, and fictitious electromotive forces in the electronic 
chiral mode fomented by the DW, on the other. An experimental realization illustrating this physics 
is proposed based on a ferromagnetic strip, which cuts the topological insulator surface into two 
gapless regions. In the presence of a ferromagnetic DW, a chiral mode transverse to the magnetic 
strip acts as a dissipative interconnect, which is itself a dynamic object that controls (and, inversely, 
responds to) the magnetization dynamics. 
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Following theoretical predictions [I] and experimental 
realizations [2 of three-dimensional topological insula- 
tors (TI's), vigorous ongoing activities in this burgeon- 
ing field are aimed at introducing spontaneous symme- 
try breaking mechanisms into the system. This could 
be accomplished by bulk or surface doping to induce 
magnetism or superconductivity in the parent (essen- 
tially free-electron) TI, or by a heterostructure design 
wherein symmetry breaking is instilled at the TI surface 
by a quantum proximity effect. We are following the lat- 
ter route, considering an insulating ferromagnetic layer 
(MI) capping the bulk TI, such that the TI surface states 
are exchange coupled to the collective magnetic moment 
of the MI. Previous theoretical investigations of a simi- 
lar TI/MI heterostructure were concerned with current- 
induced spin torques experienced by a monodomain MI 
[I], Gilbert damping by a doped TI [3], electric charging 
of magnetic textures [I], and the rectification of charge 
pumping by a monodomain precession [3], all in case of 
a well-defined spatially uniform sign of the time-reversal 
symmetry breaking gap in the TI. The essential physical 
ingredient underlying the key ideas in these papers is the 
axion electrodynamics [3] associated with the TI [5] , with 
a quantized magnetoelectric coupling that is odd under 
time reversal. In this Letter, we are interested in salient 
features associated with dynamic magnetic textures that 
imprint a spatially inhomogeneous gap onto the TI sur- 
face states, both in regard to its magnitude and sign. The 
latter, in particular, engenders electronic chiral modes at 
the magnetic domain boundaries |T] , whose hydrodynam- 
ics become intricately coupled with magnetic precession. 

According to the spin-charge helicity of the TI elec- 
tronic states, the spin-transfer torques acting on the MI 
are locked with the self-consistent electronic charge cur- 
rents in the TI. These currents, in turn, can respond to a 
combination of electromagnetic fields and fictitious forces 
induced by MI dynamics, having several distinct contri- 
butions: (i) two-dimensional (2D) surface currents re- 



lated to the half-quantized anomalous Hall effect, whose 
sign depends on the orientation of the capping magnetic 
domain, (ii) persistent currents governed by the magne- 
tization texture in the capping MI layer, and (iii) Fermi- 
level chiral currents along the domain walls (DW's) sepa- 
rating regions with an opposite Hall conductance. As an 
illustrative example, we will describe how the DW po- 
sition and an internal coordinate that parametrizes its 
Bloch-to-Neel transformation are responding to a chiral 
TI current flowing along the DW. Considering the in- 
verse charge current pumped by the DW dynamics, we 
highlight a peculiar structure of the Onsager reciprocity, 
which reverses the DW magnetization as well as the chi- 
rality of the associated electronic mode. 

Our focus will be centered on a ferromagnetic DW sep- 
arating regions with out-of-plane magnetization direction 
deep into the respective domains (which is true for suffi- 
ciently thin films, e.g., CoFeB alloys [ID]). See Fig.[T]for a 
schematic. Let us treat the DW as a stiff solitonic quasi- 
1D object, parallel to the y axis, whose translational mo- 
tion and soft internal dynamics can be described by gen- 
eralized coordinates [11]. To be specific, we start with 
the following generic free energy for an isolated magnetic 
film with magnetic spin texture m(r) (|m| = 1): 

F [m] = i J d\ {A [{d x mf + (d y mf] - Km 2 z ) , (1) 

where A is the exchange stiffness parameter and K > is 
the out-of-plane anisotropy constant. A one-dimensional 
DW running along the y axis and separating magnetic 
domains with m z = ±1 at x — > ±oo, which minimizes 
free energy ([!]), is then given in polar angles by 

0(r) = 2tan- 1 e-^»)/ Ad », 0(r) = dw , (2) 

which parametrize position r = (x, y)-dependent mag- 
netization direction m = (sin# cos <f>, sin# sin <f>, cos 9). 
A dw = y/A/K is the DW width. The DW energy is 
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FIG. 1. Schematic of a domain wall (DW) in a ferromagnetic 
strip with an out-of-plane easy (z) axis anisotropy, deposited 
on the surface of a topological insulator (TI). The DW (of 
width Adw) is parametrized, according to Eq. j2j, by two soft 
dynamic coordinates: its position Xdw(t) an d azimuthal angle 
(f>dw(t). At the DW position, a;dw, the magnetization m lies 
fully in the xy plane (forming angle </>d w with the x axis). A 
chiral electron mode (of width £ <C Ad w ) formed in the TI 
under the DW carries transport current Tdw at its exit point, 
which is governed by the voltage V y applied to the TI surface 
at its entrance and the fictitious electromotive force generated 
by the DW dynamics along its length. An Onsager-reciprocal 
spin torque affects DW dynamics in the presence of 7dw 

degenerate with respect to the position Xdw and the az- 
imuthal angle 4>dw In particular, if <f> — or tt we have 
a Neel wall and if <f> = ±ir/2 a Bloch wall. In practice, 
however, the degeneracy with respect to Xdw is lifted by 
spatial pinning fields, while the degeneracy with respect 
to ^>dw by a homogeneous applied field or spin anisotropy 
(induced, e.g., by spin-orbit interactions) in the xy plane, 
which we will take account of below. For sufficiently gen- 
tle perturbations of this kind, the zero modes associated 
with Xdv,{t) and 4>dw(t) are thus converted into soft col- 
lective excitations, which are at the core of our analysis. 

We now set out to develop a self-consistent hydro- 
dynamic theory for a DW bound with a gapless chiral 
mode, which interacts with regions of incompressible TI 
Hall fluids flanking it on the sides. The introductory 
material on magnetoelectric properties of the TI, its ex- 
change coupling to the MI, and the emergence of a chiral 
mode bound to a DW is relegated to the Supplementary 
Text (ST). From the electronic-structure point of view, 
the chiral mode patches two quantum Hall regions whose 
TKNN invariant [H] changes by unity, between the val- 
ues of ±1/2 imprinted by the magnetic domains. The 
underlying magnetoelectric effect is fundamentally dis- 
tinct from the one discussed in Ref. [4], as the DW here 
coexists with the parity-anomaly point m z — 0. 

A complete picture of our coupled magneto- 
hydrodynamic system requires us to consider also the 
spin-transfer torque that is reciprocal to the DW-driven 
electromotive forces [13j . Such torque acting on the mag- 
netization is given, due to the MI/TI exchange (see ST 



for details) 

H' = J(m x a x + rriy&y) + J±_m z a z (3) 

by (within the Landau-Lifshitz phenomenology |14| ) 

S d t m\ T = (8 m H') x m = (Jcr xy + Jj_a z z) x m. (4) 

Here, cr = (cr x , o~ y , <j z ) = {(T xy ,a z ) is the TI surface spin 
density (defined by cr = (cr) , in terms of Pauli ma- 
trices cr) and S is the saturation spin density of the 
ferromagnet. It follows from the Dirac Hamiltonian, 
Hq = v(p — eA) • z x cr + eip, furthermore, that the 
in-plane spin density cr^j, is essentially equivalent to the 
charge current density, since 

j = -(S A H) = evz x cr (T xy = (eu) _1 j x z . (5) 

Eq. ([5]) is an exact identity between the total (i.e., equi- 
librium plus nonequilibrium) current density and the in- 
plane spin density of the TI electrons, which is unspoiled 
by electron-electron interactions and ferromagnetic prox- 
imity. In particular, the chiral states, which propagate in 
the y direction and have spin quantized along the x axis, 
carry a 2D equilibrium current density that is estimated 
as (see ST for details) 

j ~ (ev/( > ){8 Aw /2nhv) z x x - (e/2irh)J±d x m z y , (6) 

where we put <5dw ~ £,J±d x m z for the chiral- mode band- 
width in terms of the characteristic (spatial) chiral-mode 
width £. From a purely phenomcnological perspective, 
on the other hand, an equilibrium charge current asso- 
ciated with a smooth static texture is given, to the first 
order in general magnetic inhomogeneities, by 

j = i]J±z x Vm z , (7) 

which should be valid both near and away from the DW, 
as long as the current is analytic in the magnetic texture 
m(r). This current is time-reversal odd, mirror symmet- 
ric (in the xy plane), and divergenceless. We, further- 
more, remark that it does not contradict the well-known 
result for the electromagnetic response of Dirac electrons 
[2], which is exact only for a strictly homogeneous sys- 
tem. The phenomcnological coefficient 77 can in general 
be a function of (J±m z ) 2 , which we take to be constant 
in the limit of weak exchange J±. Comparing Eqs. ^ 
and 0, we conclude that 77 ~ e/2ith in our model, which 
is suggestive of a universal result (as long as £ <C Adw)- 

In the presence of an equilibrium texture-induced cur- 
rent, the spin torque is given by 

S d t m\ T = 5 m F T xm, (8) 

in terms of the TI free-energy functional -FV[m] engen- 
dered by the MI/TI exchange. According to Eq. Q, this 
free energy F T can be explicitly found by integrating 

5 mxy F T = Jcr xy = (J/eu)j x z = (r]JJj_/ev)Vm z (9) 
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over m xy , at a fixed m z [IB] : 

F T = {riJJ^/ev) J d 2 r m • Vm z + i^[m*] , (10) 

where F' r is a functional of m z only [which therefore 
must derive entirely from the J± exchange in Eq. (p3j)]. 
To the leading order in the MI/TI exchange coupling 
J_i_ , F' T contributes merely to the out-of-plane anisotropy 
K in Eq. 0, which can be absorbed by a redefinition 
K — » T^f + C( Jjj = if*. Higher-order terms in F' T , in- 
cluding those that depend on spatial inhomogeneities in 
m z , would appear only at order jj_ (while cubic terms 
are prohibited by the time-reversal invariance). The 
leading-order MI/TI exchange coupling thus produces an 
anisotropy cx JJ±, which enhances tendency to form 
magnetic textures (such as skyrmion lattices), and a 
texture-independent (easy-axis) out-of-plane anisotropy 
oc J^_, corresponding to the first and second terms in 



Eq. (10), respectively. 

In addition to the equilibrium current density Q, 
there are also surface currents driven by the real and ficti- 
tious electromagnetic fields and the current carried by the 
gapless chiral mode. The latter may result in dissipation 
if connected to reservoirs (such as ungapped TI regions) . 
All these currents contribute to the torque Q. [If 7 dw is 
the ID chiral current, the corresponding 2D current den- 
sity in the y direction is jd w ~ IdwS(x — Xd w ), which is 
localized on the scale of the chiral-mode width £.] In par- 
ticular, as discussed in the ST, the torques arising from 
the effective electric-field-induced currents correspond to 
the Chern-Simons action associated with the effective 3- 



potential An = An 



[with Afj, = (ip, 



-A) denoting 
(J/ev)m x z, 



the physical and a M = (0, —a), where a 
the exchange-induced contributions]. 

Henceforth focusing on the configuration sketched in 
Fig. [T] a finite-length DW cuts across a ferromagnetic 
strip connecting semi-infinite gapless 2D reservoirs flank- 
ing its sides. In this case, the reservoirs provide an equili- 
bration and dissipation mechanism for the dc transport. 
In particular, at low frequencies, the chiral current is 
given by the Landaucr-Biittikcr formula [17] 



Idw = 9q I dy£ y = g Q V v 
Jdw 



(11) 



for the current in the y direction in response to the to- 
tal effective field £ y applied along the DW length (the 
magnetic strip width) . V y = V y +v y is the corresponding 
effective voltage (V y applied and v y induced by magnetic 
dynamics) and <?q = e 2 /h is the conductance quantum. 
The current idw in Eq. |TT] ) is defined at the exit point of 
the chiral mode and, concerning the applied voltage V y , 
only the effective electric field along the DW wire and the 
chemical potential applied to the entrance point of the 
chiral mode need to be included. The chemical potential 
at the exit point of the chiral mode, on the other hand, 




FIG. 2. The Onsager reciprocity relates voltage-induced DW 
dynamics (via spin torques) in the top panel [Eq. ( |14[ )] to 
the magnetization-dynamics-generated current (via fictitious 
electromotive force) in the bottom panel [Eq. \16\ ]. Note 
that the DW's in the bottom panel are mapped back onto 
their time-reversed parents in the top panel by a n rotation 
in the xy plane. This means that Q pumped by </>d w for the 
right chiral mode is the same in both panels. The left DW is 
treated as pinned (and thus magnetically inert) in our treat- 
ment. However, when the electron-electron interactions are 
taken into account, electrostatic charge imbalance produced 
by fictitious forces near one DW could induce currents also 
along the other DW, making such double-DW system gener- 
ally coupled. 



has no effect on the current (at both the exit and the 
entrance) in the corresponding DW. We emphasize that 
the current entering the chiral mode can generally be dis- 
tinct from Idw In particular, the dynamically- induced 
voltage v y as well as the voltage due to an electric field 
applied along the DW do not affect the entrance current, 
which is fully governed by the chemical potential applied 
at the respective lead. In this case, any imbalance be- 
tween the currents at the ends of the DW is absorbed by 
the gapped 2D regions flanking the chiral mode [in ac- 
cordance with the effective magnetic field B z = B z + b z 
[2], where b z = z • V x a = — (J/ei>)V • m is the texture- 
induced field], which we schematically sketch in the bot- 
tom panel of Fig. [2] Since the currents entering and 
exiting each individual DW thus depend very sensitively 
on the electrostatic considerations concerning the break- 
down of the effective electrochemical potentials into the 
electric and chemical counterparts, we will focus on the 
noninteracting (i.e., well-screened) electrons driven by a 
combination of a c/iemica/-potential bias at the leads and 
magnetization dynamics along the DW. 

We are now fully equipped to derive the equations of 
motion for the collective soft DW coordinates Xd-w(t) and 
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0dw(i) that parametrize the DW position and internal 
structure according to Eq. pi). In the presence of the 
(equilibrium and nonequilibrium) current-induced spin 
torques [corresponding to Eqs. (ITo| and (ITT]), respec- 



tively] as well as a uniform magnetic field h applied in 
the y direction, the full Landau-Lifshitz-Gilbert (LLG) 
equation |14j for the magnetization dynamics becomes 



5(1 + amx)9 t m = m x H* 



(12) 



where a is the dimensionless Gilbert damping constant 
and the total effective field [including the usual Larmor 
piece H c ff = —S m Fo and the spin torques] is given by 

H* =Adlra + K*m z z + hy 

+ 77* (zd x m x - yid x m z ) - j*n8{x - x dw ) ■ (13) 



Here, 



77* 



r]JJ±/ev (with r\ ~ e/2-KK), according to 
Eq. pl)| ), and J* = (J/ev)I dw , according to Eq. I dw 



is the average transport current flowing under the DW 
along the y axis 18 . 77* and j* thus parametrize the equi- 
librium and nonequilibrium spin torques, respectively. 

The equations of motion for the generalized coordi- 
nates {qt} = {idwi^dw} are derived from Eqs. |i"2"[ ) and 
(13) by integrating/ d 2 r<9 gi m-(mx [Eq. (12)]) fllj. upon 
substitution of ansatz ([2]). The key underlying phys- 
ical assumption in this procedure is that the internal 
DW structure is dominated by the A and K* terms in 
Eq. (13), such that it has a fixed width Adw ~ y/ A/K*, 



while the dynamics of slow variables qi are governed by 



the other terms in Eq. (13). Carrying out this program, 



we get (after a somewhat tedious but straightforward cal- 
culation) the following simple equations: 



x dv 



Here, 



J. sm</) dw • a . 

»dw = Xdw ■ (14) 



{4 + a 2 )S 



2A 



dw 



U = -£ d 4d*, F = ^ SiU ^ dw + ^^dwTT COS 0dw (15) 

is the generalized force (per unit of DW length L) ther- 
modynamically conjugate to the angle 4> dw . Since the 
domain wall is not pinned in the x direction, the force 
— <9 Xdw F conjugate to x dw vanishes in our model. The en- 
ergy dissipation P = -{d^ dv F)(j) dw - (d Xdv/ F)x dw associ- 
ated with magnetic dynamics (in the absence of trans- 
port current j») is thus guaranteed to be positive in 
an out-of-equilibrium situation when a > 0. The spin- 
torque-driven DW dynamics in Eq. ( 14 ) reminds us of 
a dc Josephson effect (Q cx simp). It is, in particular, 
noteworthy that the equilibrium and nonequilibrium spin 
torques add up, such that the latter can be formally ab- 
sorbed into a redefinition of 77*: 77* — > 77* + 2j«/ir. In the 
absence of the applied field, h = 0, the dynamics would 
thus settle down at <f) dw = or tt (a Neel wall) , for 77, 
(corresponding to the ordinary or n Josephson junction, 



respectively) . In the absence of spin torques but a finite 
field h along the y axis, the dynamics (that are over- 
damped as 4>dw oc fcf,) would flow towards (/> dw = ±7r/2 
(a Bloch wall), for ft. ^ 0, which corresponds to the lowest 
magnetostatic energy. 

Supplementing Eqs. ( 14 ) with the Onsager reciprocity 
principle [Tj5], dictates how the DW dynamics induce 
transport current along the chiral mode. (See Fig. [2]) 
To infer this, consider a voltage V^-induced current: 
3* = (dQJ/ ev )Vy From Eqs. (14 1 and |T5| ), which de- 
scribe how this voltage induces dynamics (x dw , <^dw)> we 
recover their Onsager (time-reversed) counterpart in the 
charge sector: 



Q 



9QJ 



a sin 0dw 



ev 2X dw {4 + a 2 )S 



9qJ_ 

ev 



Lcb dw sin cj) dv 



_ 9qJ 



ev 



Ld t m x {x dw ), (16) 



where on the second line we dropped the term that is 
diagonal in the charge sector and thus outside of the re- 
ciprocal reasoning [20] • In the final equality of Eq. (16), 



we recognize exactly the Landauer-Buttiker formula ( 11 ) 



for the magnetization-dynamics-driven charge current. It 
is crucial to notice that the DW chirality flips under time 
reversal, as illustrated in Fig. [2] The charge Q in Eq. ( [To]) 
pumped by the DW dynamics in the top panel of Fig. [2] 
thus enters the reservoir that is opposite to the one where 
the voltage V y is applied, as must be since certainly 
induces the current only downstream of the chiral mode. 
This proves internal consistency of our theory. 

In summary, we developed a self-consistent hydrody- 
namic description of a magnetic DW bound with its 
parity-anomaly chiral electron mode. DW dynamics 
parametrized by slow variables x dw and 4> dw share sim- 
ilarities with ac/dc Josephson relations for charge and 
phase, respectively. In particular, the DW switches be- 
tween two types of Neel walls (corresponding to and 
7r junctions) depending on the sign of the spin torque, 
and two types of Bloch walls depending on the sign of 
the applied field. Reciprocally, the chiral transport is 
pumped by the DW dynamics, in accord with fictitious 
gauge fields along the DW length. This coupled system 
provides a ballistic electron interconnect, which can be 
imprinted onto TI surfaces and dynamically controlled by 
magnetic fields, opening rich possibilities for "magnetic 
lithography" of electronic nanostructures on TI surfaces. 
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ELECTROMAGNETIC RESPONSE OF TOPOLOGICAL INSULATORS 

A simple strong topological insulator surface has a single Dirac cone of gapless electronic excitations pp. In the 
presence of a mirror-symmetry plane normal to the surface, and neglecting any warping effects, the effective low-energy 
Hamiltonian is then given by H Q — vp ■ z x <r + Aa z . Here, p = (p x ,p y ) — — ifiV is the 2D canonical momentum in 
the plane of the surface, z normal unit vector, <x = (a x , a yi a z ) vector of Pauli matrices for effective electron spin, v 
Fermi velocity, and A (parity-breaking) mass term due to a time-reversal symmetry-breaking mechanism, which we 
take to be infinitesimal for the moment. The linear Dirac dispersion is valid up to a momentum cutoff p c ~ E g /v, 
where E g is the bulk spin-orbit gap of the TI. Minimally coupling electrons to the gauge 3-potential = (ip, — A), 

H — > v(p — eA) • z x <x + Aa z + eip , (1) 

we recall the long- wavelength low- frequency charge response [2]: 

f = sgn(A)(e 2 /4h)e^F af! . (2) 

Here, e' 1 "' 9 is the antisymmetric Levi-Civita tensor with e 012 = 1, F a p = d a Ap — dp A a is the 3x3 field strength 
corresponding to the in-plane electric field E = —d t A — V</? and the normal magnetic field B z = z • V x A. 
j = a is the 2D surface charge density and j ' are respectively the x, y current-density components. The continuity 
equation d^j^ = is, according to Eq. guaranteed by the Faraday's law of induction (the Bianchi identity): 
z • (<9tB + V x E) = 0. We can see this explicitly after rewriting Eq. ^ in terms of the electromagnetic fields: 

a = -~sgn(A)g H B z , j = sgn(A)g H z x E , (3) 

where gn =■ e 2 /2h = i?q/2 is the Hall conductance given by half of the conductance quantum <?q = e /h. Integrating 
out electrons, would produce the Chern-Simons action for the electromagnetic field [2]: 

Scs = - (ffo/4) J dtd 2 rsgn(A)e^A^d a A . (4) 

One can easily verify that Eq. Q is in one-to-one correspondence with the charge response p]), as = —Sa^S. 
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The TI surface response ^ is tied up with the axion electrodynamics of the TI bulk [3J : The full time-reversal- 
invariant electromagnetic Lagrangian in the presence of a TI bulk is given by £em = £o — #E • B, where Cq is the 
ordinary (vacuum) Maxwell's Lagrangian and 6*E • B is the topological axion term with 8 = sgq/2, where s — sgn(A) 
inside the TI and outside. The Maxwell's equations are then modified as follows: V • E = p + V6* • B and 
VxB = <9 t E + J — (dtO + V# x E), where p and J are the ordinary external three-dimensional charge and current 
densities (entering Co). For a time-independent 6, we thus find an additional 2D charge density a = —sgqB z /2 and 
current density j = sgQZ x E/2 at the TI surface, in accord with Eqs. Finally it is useful to note that the sign s 
ambiguity of 9 inside of the TI when A — > is relieved by the invariance of the full Lagrangian upon 9^-9 + ngQ , 
where n is an arbitrary integer [3J. 



MAGNETIC-INSULATOR/TOPOLOGICAL-INSULATOR EXCHANGE INTERACTION 

Coupling the TI surface to a MI film with a dynamic magnetic (spin-density) texture S(r, t) — Sm(r, t) will induce 
time-reversal symmetry-breaking terms in the effective Hamiltonian H = Hq + H', with the axially-symmetric local 
coupling given by the anisotropic exchange: 

H = J(m x a x + niydy) + J±m z a z . (5) 

Note that there is no a priori reason to assume that J = Jj_ as in Refs. based on symmetry considerations. The 
in-plane exchange induces an effective vector potential while the out-of-plane exchange opens a gap: 

a = (J/ev)m x z and A = J±m z , (6) 

respectively. It thus follows that a dynamic magnetic texture induces a charge response given (away from the m z = 
regions) by [3] 

a = sgn(A)^V-m, j = -sgn(A)^zx9 i mxz, (7) 
lev lev 

where V should be understood as the 2D differential operator (d x ,d y ) tangential to the TI surface. Note that there 
is a sign ambiguity of the Hall response Q when m z = (and thus A changes sign), a point of central interest in 
this work. 



CHIRAL ELECTRON MODE 



The TI is gapped away from the DW, while the region under the DW supports a chiral electron mode. It can be 
physically revealed by applying a uniform electric field E along the y axis (for a configuration sketched in Fig. 1 of 
the Main Text), which induces a Hall current j x = —sgn(m z )gQE/2 (supposing the magnetization is kept static and, 
to be specific, taking J± > here). This, in turn, induces a nonequilibrium ID electron charge density A along the 
DW, which is found by the continuity equation: 

<9*A = J dxdtp = j x - - j x ,+ = gQE , (8) 

where ± designate domains at x ^ £dw, respectively. Here, we have invoked the Hall response ^ away from 
the DW (which would break down in close vicinity of the point m z = 0). On the other hand, recall that a ID 
chiral mode is excited by an electric field E as follows: the acceleration p = eE induces a charge density increase 
d t \ = ep/2irh — gqE, in agreement with Eq. ([8l). 

In order to develop a microscopic understanding of the chiral mode and evaluate the relevant energy and length 
scales, we can explicitly solve the eigenvalue problem for the Hamiltonian 

H = -ihv(a x d y - a y d x ) + C,xa z , (9) 

near the DW region where the gap A(x) = (x (£ ~ Ji/Adw) vanishes. The problem is fully analogous to calculating 
edge states of the integer quantum Hall effect [5]- As the Hamiltonian is translationally invariant in the y direction, 
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we look for energy e eigenvalues of the form ip(x, y) 
can be easily found as [6] 



<p(x), where the x-dependent part <p for the chiral mode 



<p ± = e-K/aftrO*.** n ± , e 



±VPy 



(10) 



u± here are the spinor eigenstates of a x with eigenvalues ±1, respectively. Assuming, for concreteness, that £/u > 0, 
only the oc u + eigenstate with energy e = vp y is physical, being localized as a Gaussian oc e -(C/2fo>)* near the origin. 
The characteristic width of the chiral mode is thus £ = \JtwjQ: the sharper the DW the narrower the chiral mode. 
The corresponding gap between the chiral zero mode and the dispersive gapped states from the adjacent 2D regions 
can be estimated as 



(11) 

The width of the chiral mode £ corresponds to a narrow region near Xdw , which we assume to be much smaller than 
the DW width Ad w in the ferromagnetic film (typically 10-100 nm [7]). Accordingly, in addition to the 2D charge 
density captured by Eq. (|7|), the chiral mode contributes ID charge density 



dtHy) = 9q£ v = 9QY ■ 



E 



J 



z x d t m(x dw ,y) 



(12) 



For a DW parametrized by Eqs. (2) of the Main Text, dtm(xdw, y) oc z, in the case of a sliding motion with x dw 7^ 
and </>dw = 0, such that the DW dynamics-induced charge density ( 12 1 vanishes. If 0dw ^ 0, on the other hand, a finite 



electromotive force (EMF) oc dtm x is generally generated by the DW dynamics. In order to account, furthermore, 
for the ID charge hydrodynamics along the DW, we would have to substitute d t on the LHS of Eq. (12 1 with the 
advective derivative: dt —> D t 



dt +v*d y . here is the velocity associated with the chiral-mode dispersion e = v*p y : 



u* = v(l + ^/), where gf parametrizes the (Luttinger- liquid) strength of the electron-electron forward scattering [5]. 
Since the current carried by an excess density A is idw = ^*A, Eq. (12) gives 



dt 



E y + —d t m x (x dw ,t) 
ev 



(13) 



neglecting the aforementioned advective corrections oc v*d y . The divergence of a dc current response to electric field 
E y reflects the absence of dissipation for a chiral mode, which is essentially insensitive to disorder. This changes 
profoundly, if we connect the chiral mode to Fermi-liquid reservoirs, as discussed in the Main Text. 
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